Conservation genetics of the yellow-bellied toad (Bombina variegata) and the common lizard (Zootoca vivipara) in the Italian Alps

Global change is heavily affecting Alpine ecosystems in term of both climate warming and anthropization and its effects have been already demonstrated for many different taxa. However, understanding the genetic consequences on wild species caused by environmental modifications is complicated. In this thesis, I analyzed the genetic variation pattern in two vertebrate species whose distribution and persistence across the Italian Alps could be, or already have been, affected by changing climatic conditions and human pressures for assessing their conservation status. I collected semi-invasive samples of the yellow-bellied toad Bombina variegata and of the common lizard Zootoca vivipara in the central-eastern part of Alpine chain for performing molecular analyses and subsequent statistical inferences. Different types of genetic data were used for different aims, such as mitochondrial and nuclear sequences for confirming the ESU status of a threatened lineage, microsatellite markers for evaluating genetic variability and demographic histories of wild populations or genomic SNPs for studying a major evolutionary phenotypic transition. The analyses suggested that some of the studied populations of both species suffer from reduced genetic variability and low effective population size, even if this pattern is not directly ascribable to recent anthopogenic and climatic changes. In the light of these results, however, specific conservation measure should be evaluated for these species, which are considered of least concern by the IUCN, in particular considering the predicted increase of temperature and expected modifications for their most suitable habitats

A fossil-calibrated phylogenomic analysis of Daphnia and the Daphniidae

In the post-genomic era, much of phylogenetic analyses still relies on mitochondrial DNA, either alone or in combination with few nuclear genes. Although this approach often makes it possible to construct well-supported trees, it is limited because mtDNA describes the history of a single locus, and nuclear phylogenies based on a few loci may be biased, leading to inaccurate tree topologies and biased estimations of species divergence time. In this study, we perform a phylogenomic analysis of the Daphniidae family (Crustacea: Branchiopoda: Anomopoda) including some of the most frequently studied model organisms (Daphnia magna and D. pulex) whose phylogenetic relationships have been based primarily on an assessment of a few mtDNA genes. Using high-throughput sequencing, we were able to assemble 38 whole mitochondrial genomes and draft nuclear genomes for 18 species, including at least one species for each known genus of the family Daphniidae. Here we present phylogenies based on 636 nuclear single-copy genes shared among all sampled taxa and based on whole mtDNA genomes. The phylogenies we obtained were highly supported and showed some discrepancies between nuclear and mtDNA based trees at deeper nodes. We also identified a new candidate sister lineage of Daphnia magna. Our time-calibrated genomic trees, which we constructed using both fossil records and substitution rates, yielded very different estimates of branching event times compared to those based on mtDNA. By providing multi-locus, fossil-calibrated trees of the Daphniidae, our study contributes to an improved phylogenetic framework for ecological and evolutionary studies that use water fleas as a model system.Peer reviewe

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the "Gymnomera"). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]Peer reviewe

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the "Gymnomera"). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]Peer reviewe

Phylogeography of the snake pipefish, <i>Entelurus aequoreus</i> (Family: Syngnathidae) in the northeastern Atlantic Ocean

The snake pipefish, Entelurus aequoreus, is a widespread marine species occurring in pelagic and coastal environments in the northeastern Atlantic Ocean. Recently, the snake pipefish underwent a short-lived, yet substantial, increase in abundance and range expansion into arctic waters. However, little is known about the species' population structure or if different ecotypes contributed to this outbreak. Specimens (n = 178) were sampled from 25 locations from six regions spanning 1.9 million km2. A fragment of the mitochondrial cytochrome b gene and control region was used to assess population structure and genetic diversity. Both loci showed high haplotype diversity (Hd) and low nucleotide diversity (π) over all sampled locations. A genetic signature of population expansion was evident through mismatch distributions and tests for recent population expansion (Fu's Fs, Tajima's D and R2). Effective population size analyses (Bayesian skyline plot) suggest an expansion 50-100 thousand years before present. However, we found neither significant population differentiation (analysis of molecular variance) among regions nor evidence of genetically distinct ecotypes. This lack of structure is probably due to a pelagic life style, fast development and long distance dispersal aided by ocean currents. Our work highlights the need for further research to better understand the recent outbreak and how this species may respond to future environmental challenge

Survival and divergence in a small group: The extraordinary genomic history of the endangered Apennine brown bear stragglers

About 100 km east of Rome, in the central Apennine Mountains, a critically endangered population of ∼50 brown bears live in complete isolation. Mating outside this population is prevented by several 100 km of bear-free territories. We exploited this natural experiment to better understand the gene and genomic consequences of surviving at extremely small population size. We found that brown bear populations in Europe lost connectivity since Neolithic times, when farming communities expanded and forest burning was used for land clearance. In central Italy, this resulted in a 40-fold population decline. The overall genomic impact of this decline included the complete loss of variation in the mitochondrial genome and along long stretches of the nuclear genome. Several private and deleterious amino acid changes were fixed by random drift; predicted effects include energy deficit, muscle weakness, anomalies in cranial and skeletal development, and reduced aggressiveness. Despite this extreme loss of diversity, Apennine bear genomes show nonrandom peaks of high variation, possibly maintained by balancing selection, at genomic regions significantly enriched for genes associated with immune and olfactory systems. Challenging the paradigm of increased extinction risk in small populations, we suggest that random fixation of deleterious alleles (i) can be an important driver of divergence in isolation, (ii) can be tolerated when balancing selection prevents random loss of variation at important genes, and (iii) is followed by or results directly in favorable behavioral changes

Survival and divergence in a small group: The extraordinary genomic history of the endangered Apennine brown bear stragglers

About 100 km east of Rome, in the central Apennine Mountains, a critically endangered population of &sim;50 brown bears live in complete isolation. Mating outside this population is prevented by several 100 km of bear-free territories. We exploited this natural experiment to better understand the gene and genomic consequences of surviving at extremely small population size. We found that brown bear populations in Europe lost connectivity since Neolithic times, when farming communities expanded and forest burning was used for land clearance. In central Italy, this resulted in a 40-fold population decline. The overall genomic impact of this decline included the complete loss of variation in the mitochondrial genome and along long stretches of the nuclear genome. Several private and deleterious amino acid changes were fixed by random drift; predicted effects include energy deficit, muscle weakness, anomalies in cranial and skeletal development, and reduced aggressiveness. Despite this extreme loss of diversity, Apennine bear genomes show nonrandom peaks of high variation, possibly maintained by balancing selection, at genomic regions significantly enriched for genes associated with immune and olfactory systems. Challenging the paradigm of increased extinction risk in small populations, we suggest that random fixation of deleterious alleles (i) can be an important driver of divergence in isolation, (ii) can be tolerated when balancing selection prevents random loss of variation at important genes, and (iii) is followed by or results directly in favorable behavioral changes.Additional co-authors: Claudio Groff, Ladislav Paule, Leonardo Gentile, Carles Vilà, Saverio Vicario, Luigi Boitani, Ludovic Orlando, Silvia Fuselli, Cristiano Vernesi, Beth Shapiro, Paolo Ciucci, and Giorgio Bertorell

No evidence for genetic sex determination in Daphnia magna

Mechanisms of sex determination (SD) differ widely across the tree of life. In genotypic sex determination (GSD), genetic elements determine whether individuals are male or female, while in environmental sex determination (ESD), external cues control the sex of the offspring. In cyclical parthenogens, females produce mostly asexual daughters, but environmental stimuli such as crowding, temperature or photoperiod may cause them to produce sons. In aphids, sons are induced by ESD, even though GSD is present, with females carrying two X chromosomes and males only one (X0 SD system). By contrast, although ESD exists in; Daphnia; , the two sexes were suggested to be genetically identical, based on a 1972 study on; Daphnia magna; (2n=20) that used three allozyme markers. This study cannot, however, rule out an X0 system, as all three markers may be located on autosomes. Motivated by the life cycle similarities of; Daphnia; and aphids, and the absence of karyotype information for; Daphnia; males, we tested for GSD (homomorphic sex chromosomes and X0) systems in; D. magna; using a whole-genome approach by comparing males and females of three genotypes. Our results confirm the absence of haploid chromosomes or haploid genomic regions in; D. magna; males as well as the absence of sex-linked genomic regions and sex-specific single-nucleotide polymorphisms. Within the limitations of the three studied populations here and the methods used, we suggest that our results make the possibility of genetic differences among sexes in the widely used; Daphnia; model system very unlikely

Data from: Defining conservation units in a stocking-induced genetic melting pot: unravelling native and multiple exotic genetic imprints of recent and historical secondary contact in Adriatic grayling

The definition of conservation units is crucial for the sustainable management of endangered species, though particularly challenging when recent and past anthropogenic and natural gene flow might have played a role. The conservation of the European grayling, Thymallus thymallus, is particularly complex in its southern distribution area, where the Adriatic evolutionary lineage is endangered by a long history of anthropogenic disturbance, intensive stocking and potentially widespread genetic introgression. We provide mtDNA sequence and microsatellite data of 683 grayling from 30 sites of Adriatic as well as Danubian and Atlantic origin. We apply Bayesian clustering and Approximate Bayesian Computation (ABC) to detect microgeographic population structure and to infer the demographic history of the Adriatic populations, to define appropriate conservation units. Varying frequencies of indigenous genetic signatures of the Adriatic grayling were revealed, spanning from marginal genetic introgression to the collapse of native gene pools. Genetic introgression involved multiple exotic source populations of Danubian and Atlantic origin, thus evidencing the negative impact of few decades of stocking. Within the Adige River system, a contact zone of western Adriatic and eastern Danubian populations was detected, with ABC analyses suggesting a historical anthropogenic origin of eastern Adige populations, most likely founded by medieval translocations. Substantial river-specific population substructure within the Adriatic grayling Evolutionary Significant Unit points to the definition of different conservation units. We finally propose a catalog of management measures, including the legal prohibition of stocking exotic grayling and the use of molecular markers in supportive- and captive-breeding programs

nABCtemolo_4pops_SSRS_I

Microsatellite data set analysed in AB